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Abstract. The symplectic tomograms of 2D Hermite-Gauss beams are found 
and expressed in terms of the Hermite polynomials squared. It is shown that 
measurements of optical-field intensities may be used to determine the tomograms 
of electromagnetic-radiation modes. Furthermore, entropic uncertainty relations 
associated with these tomograms are found and applied to establish the 
compatibility conditions of the the field profile properties with Hermite-Gauss 
beam description. Numerical evaluations for some Hermite-Gauss modes 
illustrating the corresponding entropic uncertainty relations are finally given. 



PACS numbers: 42.50.-p, 42.50..Dv, 03.67.-a 

1. Introduction 

The paraxial optical beams propagating in media are known to obey the Schrodinger- 
like equation [TJH]. The properties of the paraxial beams described by quantumlikc 
equations were intensively discussed in both standard optics (see, for example, [2111]) 
and charged-particle-beam optics [5]. In particular, by means of such quantumlikc 
equations, the paraxial propagation of the optical-field modes in both free space and 
in media, where the refractive index is a given function of the spatial coordinates, can 
be described in terms of the Hermite-Gauss or the Laguerre-Gauss modes [5J |7] . 

Recently, the tomographic analysis of electromagnetic beams received a great deal 
of attention, since it allows for an accurate characterization of the beam mode profile 
required when direct observations are not able to determine the phase of the optical 
fields and, therefore, measure of the beam intensity in all points is, in principle, needed 
0. 

In view of a tight connection between the paraxial beam propagation and quantum 
mechanics mentioned above, it seems quite natural to employ a quantum tomographic 
approach j9] to characterize optical beam mode profiles. Remarkably, to this regard, 
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it was found that the quantum states, being the solutions to the Schrodinger equation, 
can be mapped by means of the Radon transform onto the standard probability 
distributions called symplectic [TO] and optical [pj tomograms. The tomograms were 
shown [llj to obey the constrains which correspond to entropic uncertainty relations 
[T2] for Shannon entropy [13] and Renyi entropy [T3] of quantum states. 

In this paper, we construct symplectic tomograms of two-dimensional (2D) 
optical modes and discuss their properties. Furthermore, we introduce the concept 
of tomographic entropy of such modes and show for them the existence of entropic 
uncertainty relations in close analogy to the Heiscnbcrg or Robertson-Schrodinger 
uncertainty relations. Numerical results showing the capability of entropic uncertainty 
relations to characterize the spatial content of 2D optical modes are presented. 



2. Symplectic tomography of two-dimensional beams 



Let us consider a 2D (in Cartesian x and y transverse coodinates) optical mode field 
of wavelength A which satisfies the following paraxial wave equation [TJ [2] : 

A 2 



iX 



dip{x,y,z) 



- V_l^ + U(x, y, z)ip(x, y, z), 



(1) 



where = 



.,, A = X/2tt, U oc —Sn(x,y,z) (n being the refractive index 



dz 

d 2 St 

dx 2 dy 2 

close to the z propagation direction), and ip(x,y,z) is the complex amplitude of 
electromagnetic field. In this equation, the longitudinal coordinate plays the role 
of time and A plays the role of the Planck's constant. The refractive-index profile 
provides the possibility to control the modes. In fact, the refractive index, being 
dependent on the transversal and longitudinal coordinates, determines "the potential 
energy" as a function of the position and time in the corresponding Schrodinger-like 
equation. 

Since we assume a real U, the following normalization condition can be imposed: 



\ip(x,y,z)\ 2 dxdy = 1. 



(2) 



The tomographic-probability distribution (called also tomogram) is defined by the 
Radon transform |15j generalized for two dimensions as symplectic tomogram |16j and 
given in terms of the squared modulus of the Fresnel integral ( [T7] ) 

w(Xi, fil, V\, X 2 , M2, V2, z) 
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Hereafter, to simplify the notation, we do not show explicitly the z dependence of the 
tomogram and related expressions. For instance, w{X\, /ii, v\, X%, ^2, ^2, z) is replaced 
by w(Xi , Hi, v 1, X 2 , /i2 j v 2) ■ The tomogram is nonnegative function of six real variables 
X\, fa? V\i X 2 , fi>2i and v 2 and it satisfies the normalization condition 



J w(Xi,ni,vi,X 2 ,fi 2 ,v 2 )dXi dX 2 = l. 



(4) 



Tomogram |3]) can be interpreted as the probability density of two random variables 
Xi = fii x + v \Vx an( l X 2 = \i 2 y + ^2Pj/ 5 where x and y are coordinates of an 
intersection point in the transversal plane of the light ray and p x and p y are small 
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angles determining the unit direction vector parallel to the light ray propagating along 
the fiber axis. 

The tomographic-probability density determines the modulus and phase factor of 
the mode function tp(x, y) due to the inverse relation 

1 



4tt 2 



w(x 1 ,fi, 1 ,v 1 ,X2,(j,2,v2)$(.vi -x + x')5(v 2 - y + y') 



x exp 



- (2X X - fi^x + x x ) + 2x 2 (y + y')y 2 - ^ 2 {y + y')) 



dXi dX 2 dfii dfj, 2 dv\ dv 2 . 

(5) 



Tomogram (|3| has the homogeneity property 

w(AiXl, Ai/xi, Ai^i, \ 2 X 2 . X 2 fi 2 , \ 2 v 2 ) = 1 w(X 1 ,[i 1 ,v 1 ,X 2 ,n 2 ,i' 2 ), (6) 
and it can be expressed in terms of optical tomogram depending on four real variables 
w op t(Xi 7 6\, X 2l 9 2 ) 



4-7T 2 1 sin 6i sin 9 2 1 



ip(x,y) exp 



j (cot + XI) + cot 9 2 (if + XI) - - 

2 V sin 0i sin #2 



dx dy 



In fact, due to the definitions ([3| and one has 

w op t(Xi,0i,X 2 ,9 2 ) = w(Xi,cos9i,sm9 1 ,X 2 ,cos9 2 ,sin9 2 ) 

and, due to the homogeneity property |6l, 

1 
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Another tomogram called Fresnel tomogram of the light mode in the optical fiber is 
given by the integral transform 18J 

2 



w F (X 1 ,v 1 ,X 2 ,v 2 ) = 



1 1 
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ip(x,y) exp 



iix.-x) 2 , i{x 2 - y y 
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2v 2 
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(10) 

Fresnel tomogram wp can be obtained from symplectic tomogram by putting fj,i — 1 
and fi 2 = 1, i.e., wf(Xi, v\, X 2 , v 2 ) = w(Xi,l,vi,X 2 ,l,i> 2 ). Furthermore, the 
following relationship holds: 



w(X 1 ,fjL 1 ,V 1 ,X2,fl2,V 2 ) 



1 

\f*lH2\ 



Xl 2/1 X 2 v 2 
Ml ' Ml ' M2 M2 



(11) 



In view of mutual tomogram relations ([9]) and (10 1, one can find the mode function 
il>(x,y) using ^ and express this function either in terms of optical tomogram or in 
terms of Fresnel tomogram. 
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3. Tomographic entropies of optical beams 

In this section, we consider in detail tomographic entropies associated with optical 
beams following the general formalism developed in [13]. There exists the Shannon 
construction [13] of entropy associated to the probability distribution function P(n) 
of a discrete variable n 

H = -^P(n)lnP(n), (12) 

n 

which can be easily extended to the classical probability distribution function of a 
continuous variable. In particular, we can use the tomographic-probability distribution 
(J3j to introduce the following tomographic entropy for an optical beam: 

H(ni,vi,H2,v 2 ) = - J w(Xi, (ii, vi, X 2 , ti 2 , V2)tnw(X 1 , fi 1 , vi, X 2 , H2, V2) dXi dX 2 . 

(13) 

The above entropy is a new information characteristic of the optical-beam profile in 
media, which yields also the optical tomographic entropy of the optical beam 

H opt (9 1 ,9 2 ) = - [ w opt (X 1 ,9 1 ,X 2 ,9 2 )]nw opt (X 1 ,9 1 ,X 2 ,9 2 )dX 1 dX 2 (14) 

and the Frcsnel tomographic entropy 

H F (vi,v 2 ) = - J w F (Xi,vi,X 2 ,v 2 )lnw F (Xi,is 1 ,X 2 ,v 2 )dXidX 2 . (15) 

All three entropies ([13])— ( 15 > are mutually related. 

An interesting property of the tomographic entropies is connected with the 
uncertainty relation associated to the light-beam intensities related to the mode 
function ijj(x,y) and its Fourier transform 

i>(p x , Py ) = (27T)" 1 J ^(x,y)exp-^ x+ P^ dxdy. (16) 

This entropic uncertainty relation reads (see, e.g., |12j ) 

\^(x, y)\ 2 In \^(x, y)\ 2 dx dy- [ \$(x, y)\ 2 In \${x, y)\ 2 dx dy > 21n(7re). (17) 



The entropic uncertainty relation was generalized for the symplectic tomographic 
entropy as well as for the optical and Fresnel tomographic entropies [18j . 

To consider, for example, the optical tomographic entropy, we introduce the 
function 

R{9 U 9 2 ) = H opt {9 u 9 2 ) + F opt (0i + tt/2, 2 + vr/2) - 2hu%e), (18) 



where H opt (6i, 9 2 ) is given by ([14 1. 

According to the new entropic uncertainty relations |18j . the function R{9i,9 2 ) 
must be nonnegative for all values of the angles 9\ and 9 2 , i.e., 

R(9 1 ,9 2 )>0. (19) 

This means that, if one measures the modulus and phase of the mode function ^{x, y) 



by any method, the results of the measurement yield also the function ([18| which 
must be nonnegative. Nonnegativity of this function for all the angles 9\ and 9 2 can 
serve as an extra control of accuracy of the measurements. 
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4. Tomography of 2D Hermite Gauss beams 

We want to illustrate here the tomographic properties of a two-dimensional Hermite- 
Gauss (HG) beam of order (n, m) described by a wave function VWi^i; x 2)- These 
beams are of a particular interest since they form a complete set of solutions of 
the paraxial equation in the case of a linear medium whose refractive index is 
a quadratic function of x\ and x 2 or in vacuo (note that in ([l]) we have replaced x 
and y by x\ and x 2 , respectively). Any other solution can be written, in principle, 
as an expansion in terms of HG beams over all the indices n and m. The HG 
beams play an important role in the design of spherical resonators since they represent 
the resonators' modes having the same wavefronts as a Gaussian beam but different 
amplitude distributions. 

Let us write the HG beam of the order (n, m) in the following normalized form: 



tpnm(xi,X 2 ) 



'a/2: 



•< - 2 



CO 



exp 



(20) 



where Co is the width assumed to be the same along the two transverse directions X\ 

I 2 . 

and xr> and N„ m = , / — f, — . — . , — is the normalization factor. 

V Tro-$n\m\2 n + m 

The 2D tomogram w nrn of HG beam of the order (n, m) is given by the following 
six-dimensional expression: 

1 1 
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(21) 



Substituting the HG field given by Eq. (20 1 into Eq. (21 1 and taking into account the 
integral relation 



H n (ay)exp -{y-(3)' 



dy = \/7T (l 



2W 2 tf 



af3 



Vl - a 2 



(22) 



we obtain for tomogram the following expression: 



u>„m(Xi, X 2 , Ml) M2, fi, ^2) 
where 



1 



1 



(2tt) 2 I^iaj 



7V 2 m |7(X 1; Ml , i/O | 2 |/(X 2 , M2 , ^) | 2 , (23) 



I(X k ,n k ,v k ) = V2 (1 - a 2 ) n/2 H n 



ctkPk \ exp (-X2/4^ % ) 



fe = l,2,(24) 



and the parameters a k , (3 k , and % are given by 

y/2 „ X,, 



OL k = 



o~o\/qk 
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Ik 
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2v k 



(25) 



In view of the definition given by Eq. (21 1 inserted into Eq. (24 1 for k = 1, 2, we 
can write the tomogram w nm (Xi, X 2 , Hi, Hi, vi, v 2 ) in a form which gives insight into 
its physical meaning, namely, 



Wnm{Xl,X 2 ,Hl,H2, Vi,V 2 ) 
2 
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where the following quantities are defined: 



a k (z) = a \ I + 



Pk_ _ _2tt_ 
v k Xz k 



X 



(27) 



Now we consider a special case p\jv\ = ^2/^2 = p- 

In this case, the tomogram is a function of four phase space variables only, namely, 



Xi, X2, p, and and Eq. (26 1 can be written as follows 

•2 
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pcr(z) 



\pv\ 



exp 



a z (z)J \ pa(z) 
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pja 2 (z) nl<r 2 (z 



(28) 



z given in Eq. ( 27 I 



where, according to Eq. (27 1, we have <x(z) = <to\/ 1 + ( z /zo) and z = 27r/Ap, with 



Apart the scaling factor (\p1p2D , Eq. (28 1 represents the intensity distribution 



of the HG beam of the order (n, m) at a distance z in the transverse plane of rescaled 



transverse distances X\j ' p and X%/ [A. The first part of Eq. (28 1 gives the well-known 
dependence of the beam radius o(z) on the propagation distance z — 2ir/Xp, which is 
defined through the second part of Eq. p8] ) in terms of the ratio p = fi/v of the phase 
space variables. At z — 0, the beam radius is given by cr(0) = cr . 



The wavelength A plays the role of a spatial scaling parameter in Eqs. (27 1 and 



(28 1. We point out the formal analogy between the tomographic representation of 



HG beams and the paraxial propagation law governing the evolution of the intensity 
distribution in the rescaled transverse plane of an optical HG beam of the order 
(n,m) and wavelength A. Indeed zq given by Eq. (27 1 is the well-known confocal 
parameter or the depth of the focus of the propagating beam. When the relation 
Pi/v\ = ^2/^2 = p is not satisfied and the phase space variables assume arbitrary 
values, the simple interpretation of tomogram in terms of the intensity distribution of 
2D HG beam at a fixed distance z — 2ir/Xp does not more hold. In general, we have 
two different distances which, according to Eq. (27), are given by Z\ — 2-KV\jXp\ and 
2iTU2/Xp2, respectively. 

However, it is easy to see that Eq. (27) is separable along the two transverse 



-2 



directions X\j 1 p and X2I p and that the corresponding factors represent the intensity 
distribution of two independent ID HG beams of radii cr 1 (z 1 ) and 02^2) at the 
distances zi and z 2 - 

Therefore, the tomogram of 2D HG comprises, in general, the features of 
propagation laws corresponding to the projections of the beam along two orthogonal 
directions depending on the values of the phase space variables. 

In Figs. 1-3, we present plots of the function R(9i, 62) calculated for three different 
cases of 2D Hermite Gauss modes. One can see a similar behavior of these functions 
for different mode indices and widths. Nevertheless, all the functions are nonnegative 
according to the found uncertainty relations. 



5. Conclusions and discussions 



In this paper, we have shown that the simplectic tomography of two-dimensional 
optical beams can be employed to characterize the spatial distribution of the beam 
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tr=l <r= 4 




Figure 1. Two-dimensional entropic inequality of the Hermite— Gauss beam Hoo 
(m = 0, n = 0) as a function of the two phase angles measured in radians with 
beam width a = 1 and 4. 




Figure 2. Two-dimensional entropic inequality of the Hermite— Gauss beam Hn 
(m = 1, n = 1) as a function of the two phase angles measured in radians with 
beam width a = 1 and 4. 

cr= 1 <r= 4 




Figure 3. Two-dimensional entropic inequality of the Hermite— Gauss beam Hio 
(m = 1, n = 0) as a function of the two phase angles measured in radians with 
beam width a = 1 and 4. 



mode along the propagation direction. With a proper identification of the phase 
parameters, the tomographic map reduces to measure the optical intensity distribution 
of the optical beam in a given transverse plane. However, a full tomographic 
characterization of the beam requires, in principle, an infinite set of measurements 
of the spatial distributions of the light intensity at different points in the cross section 
of the optical beam along the propagation direction, which, in practice, is impossible 
due to the bandwidth limitations of the optical detectors. Nevertheless, we have shown 
that the tomographic approach can be employed to derive a set of entropic inequalities 
in terms of the logarithmic measurement of the beam intensity. An advantage of this 
formulation is that it allows one to catch the relevant features of the propagation 
beam characteristics from a tomographic map. Here, the concept of entropy, well 
established by statistics and later extended by Shannon to information theory, plays 
the role of an integral descriptor of the tomogram of the optical mode in terms of 
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two basic optical beam parameters, the width and Rayleigh range. According to 
the tomographic-based entropic inequality, a strictly positive function of two space 
variables 11(81,62) can be associated to the optical mode. The numerical results 
obtained for the basic Hermite-Gauss modes show clearly that this entropic-based 
optical descriptor is sensitive to the beam width and, in general, tends to increase 
for higher-order modes. It is also expected that this formulation can be extended 
to optical vortex field where the presence of several defects or vortices is a manifest 
signature of high-order optical modes. 
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